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ABSTRACT 

The deprojection of the surface brightness distribution of an axisymmetric galaxy does 
not have a unique solution unless the galaxy is viewed precisely edge-on. I present 
an algorithm that finds the full range of smooth axisymmetric density distributions 
consistent with a given surface brightness distribution and inclination angle, and use 
it to investigate the effects of this non-uniqueness on the line-of-sight velocity profiles 
(VPs) of two-integral models of both real and toy discy galaxies viewed at a range of 
inclination angles. Photometrically invisible face-on discs leave very clear signatures 
in the minor-axis VPs of the models (Gauss-Hermite coefficients hi > 0.1), provided 
the disc-to-bulge ratio is greater than about 3%. I discuss the implications of these 
hitherto neglected discs for dynamical modelling. 

Key vifords: celestial mechanics, stellar dynamics - galaxies: elliptical and lenticular 
- galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 

A fundamental problem in astronomy is determining the dis- 
tribution of instrinsic shapes of elliptical galaxies. One way 
to do this is by assuming that every galaxy is a triaxial eUip- 
soid, and then to use simple photometrical and kinematical 
information to try to constrain the distribution of axis ratios 
of their outer regions (e.g., Binney 1985; Franx, lUingworth 
& de Zeeuw 1991; Statler 1994). About a third of all ellipti- 
cals, however, have "discy" isophote distortions (Bender et 
al. 1989), steep power-law central density cusps, and show 
no obvious indications of triaxiality (Faber et al. 1997; Ko- 
rmendy & Bender 1996; see also Merritt & Quinlan 1998). 
Little is known about the discs that most probably cause 
these distortions. For example, it is unclear whether discy 
ellipticals form a continuous sequence in disc-to-bulge ra- 
tio with SO galaxies (Scorza & Bender 1995). Moreover, it 
is possible that quite large discs could lurk in those power- 
law galaxies that do not have obvious isophote distortions: 
Rix & White (1990) have shown that even the large discs 
in SO galaxies are photometrically undetectable for about 
50% of orientation angles. Therefore, rather than trying to 
determine axis ratios for these power-law galaxies under the 
assumption of ellipsoidal symmetry, it is perhaps more in- 
teresting to use kinematical information to try to determine 
(or at least constrain) their full central stellar density dis- 
tributions p{R, z) under the assumption of axisymmetry. 

The ease with which discs can be hidden can be un- 
derstood using the Fourier slice theorem (Rybicki 1987): 
in Fourier space, the surface brightness distribution of an 
axisymmetric galaxy with a given inclination angle i pro- 
vides information about the density only outside a "cone 
of ignorance" of opening angle 90° — i around the galaxy's 
symmetry axis, where i = 90° for an edge-on galaxy. Ger- 
hard & Binney (1996) have explicitly constructed a family of 
disc-like "konus densities" whose Fourier transforms vanish 



outside the cone of ignorance, and therefore whose densities 
project to zero surface brightness. Thus, the photometric de- 
projection problem is not unique for non-edge-on galaxies: 
any solution p{R, z) is degenerate to the addition of disc- 
like konus densities. Further examples of konus densities are 
given by Kochanek & Rybicki (1996) and van den Bosch 
(1997). 

Konus densities are quite difficult to construct analyti- 
cally, so a numerical approach is needed to explore the full 
range of possible solutions to the photometric deprojection 
problem. Such an approach was taken by Romanowsky & 
Kochanek (1997; hereafter RK) who used a maximum pe- 
nalized likelihood algorithm, penalizing their solutions with 
"bias functions" that influenced the shape of the final so- 
lution chosen, as well as providing the smoothness criterion 
that is necessary for any discrete numerical approximation. 
They found that the uncertainties in the deprojection were 
even greater than one would expect from the simple an- 
alytical examples of konus densities above. RK also used 
Jeans equations modelling to predict the projected second- 
order velocity moments corresponding to each solution, un- 
der the assumption of a constant mass-to-light ratio and 
a two-integral distribution function. They found that the 
second-order moments were only very weakly affected by 
the non-uniqueness of the deprojection: different solutions 
gave almost identical projected second-moment profiles. 

Real kinematical measurements, however, include more 
information than just the second-order velocity moments. 
Indeed, nowadays one routinely measures the full line-of- 
sight velocity profiles (VPs) of galaxies. In this paper, I use 
simple models to investigate whether using VPs can place 
stronger constraints on the intrinsic shapes of axisymmet- 
ric galaxies. It is organized as follows. Section 2 describes 
the deprojection problem in numerical terms and gives my 
method of solution. Like RK, I use maximum penalized like- 
lihood, but with a penalty function that trades off the ele- 
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gance of their bias functions for comprehensibility and ease 
of use. I also use an interpolation scheme that is better able 
to handle the cusped density profiles present in real galaxies. 
In Section 3, I demonstrate the effects that konus densities 
have on the VPs of simple toy galaxy models. The discs in 
real galaxies are not strict konus densities, so in Section 4 I 
use models of some real edge-on discy ellipticals to investi- 
gate whether their VPs would contain any signature of the 
discs were they viewed closer to face on. Section 5 sums up 
and suggests future work. 



2 NUMERICAL SOLUTION OF THE 
PHOTOMETRIC DEPROJECTION PROBLEM 

For projected quantities, let us choose a system of co- 
ordinates {x' ,y' , z') with origin O at the centre of the galaxy, 
the Oz'-axis parallel to the line of sight to the galaxy, and 
the Ox'-axis parallel to the galaxy's projected major (mi- 
nor) axis for an oblate (prolate) model. Real galaxies are 
close to spheroidal, so assume that we are given discrete 
measurements of the projected surface brightness distribu- 
tion I{x' , y') on an elliptical polar grid {m! , 6') with n^/ xn^/ 
vertices at 



Xij = nij cos 6j 



(1) 



Pij = qrrii sm 9j , 

where m^ runs logarithmically between the radii of the inner- 
and outer-most isophotes, 9'j is spaced linearly between 
and 7r/2, and q' (1/g') is the mean projected axis ratio of the 
galaxy assuming it is oblate (prolate). Let Sij = log I{m'i, 9'j) 
and let /S.Sij be the corresponding measurement error (as- 
sumed Gaussian). 

We seek a three-dimensional luminosity density p{R, z), 
which, when viewed at a given inclination angle i, projects 
to an acceptable fit to this Sij. We expect this p{R, z) to be 
roughly spheroidal, and to have a radial profile that is locally 
well-approximated by some power law. Hence I represent 
p{R,z) as TZ(m,9) = log p(m,^) on an rim x rig grid with 
vertices at 



Rij = nii cos 9j 
Zij — qnii sin 9j , 



(2) 



where the intrinsic axis ratio q is related to q' by q'^ = 
g^ sin^ i + cos^ i, 9j runs linearly from to it/2, and m^ 
runs logarithmically from m'l to a few times m^/ . Values 
of p{m, 6) between grid points are obtained by interpolat- 
ing linearly in TZ and exponentiating. The calculation of the 
model's projected surface brightness Sij on the grid (1) is 
outlined in the Appendix. 

Given these Sij , the goodness of fit is measured by 



n / riQf 

1=1 i=i 



A5,. 



(3) 



A reasonable model will have x — ^m'l-e' i v^/rt^/Jngt . 

There will be many distributions TZij that will have a 
value of x^ in this range, but not all of them will be ac- 
ceptably smooth. In addition, one would like the ability to 
choose among the various acceptable solutions, for example 
by preferring those with a given degree of disciness or boxi- 
ness. Given these reasonable prior prejudices (let us denote 



them as /), and an observed Sij, Bayes' theorem gives the 
relative plausibility of any set of TZij as 



p(7^|5,/)«p(5|7^,/)p(7^l/), 



(4) 



where the likelihood p(5 | TZ) = exp( — ^x^) ^nd the prior 
piTZ I /) expresses our preconceptions (smoothness, isocon- 
tour shape) about the form of p{R, z). Taking the logarithm 
of equation (4) yields an expression for the penalized log- 
likelihood 



C = ~-x'+Pi[TZ], 



(5) 



where the penalty function Pi[TZ] = logp(7?- j I) penal- 
izes those solutions that do not fit our smoothness and 
isocontour-shape requirements. 



2.1 The penalty function 

One of the simplest ways of penalizing non-smooth solutions 
is by using the mean-square second derivative: 
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where Cm and Cg are the weights given to radial and angular 
smoothness and I have omitted some uninteresting constant 
factors. Sensible values for the weights can be estimated as 



follows. Consider a typical galaxy with p 



in the 



outer parts, rolling over to p ~ m^^ in the centre. Thus 
972-/9 log m changes from —4 in the outer parts to —1 near 
the centre, and so the first sum in equation (6) must be at 
least S'^/rim- Choosing Cm ~ —{2nm>ngi)^''^nm/^mj SiU in- 
crease in the RMS change in dTZ/d log m by an amount Am 
is considered as bad as an one-sigma increase in x? of 
{2nmingiY''^ . I choose Am = 3 for all the results in this 
paper. Similarly, for the angular smoothing weight I choose 
Cg — —{2nm'ngr)^'^ng/\^. Since not much is known about 
galaxies' angular profiles, it is not so clear how to choose 
\g a priori, but experiments with real galaxies show that 
Ag ~ 0.5 is about right. 

The shape (degree of disciness or boxiness) of p{m, 9) 
can be measured by its cos 49 Fourier coefficients, defined as 



C4{m) 



7r/2 



p(m, 6*) cos 46* de*. 



(7) 



-tt/2 



similar to the various definitions often used to quan- 
tify isophote shapes (e.g., Jedrzejewski 1987; Bender & 
MoUenhoff 1987). At radii where the model is locally discy 
(boxy) C4 will be positive (negative). The penalty function 



shape 



\TZ] 



Tlfji 



1/2 
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C4,{mi) — d 
Ad 



(8) 



favours solutions with C4 ~ d±Ad, provided the observations 
do not place strong constraints on C4. If on the other hand 
the observations constrain C4 to lie outside the range d ± 
Ad, Pshape biases the solution towards C4 — d, with the 
parameter Ad controlling the strength of the bias. 
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A final, not-so-obvious, necessary constraint on the 
shape of the density isocontours is that dTZ/dO be non- 
negative aX 6 = and 7r/2. I use 



C, 



Pnn = ^^^ (max(0,7^..2 -7^.,l)) + 



c, 



— ^ 2_. ( max(0, 7?.i,„(, - T^i.no-i))" 



(9) 



to penahze each radial grid point that does not satisfy these 
conditions by an amount Cnn = —A^rim'ngi . 

Taking equations (6), (8) and (9) together, the full 
penalty function is Pi = Psmooth + -Pshape + Pnn- 

2.2 Finding the best solution 

The following procedure is used to find the density p{m, 6) 
that maximizes the penalized log- likelihood C. There may 
well be more efficient schemes, but this one is easy to im- 
plement and works well. First the parameters {L, mo, a, 13) 
of the distribution 



p(m) 



L 



./5-3 



47rgB(3 - a, 3 - /3) m"{mo + m)l^' 



(10) 



that minimize x^ are found using the Levenberg-Marquardt 
algorithm (Press et al., 1992). The function B(-,-) in (10) 
is the usual Beta function. Then the following simulated 
annealing scheme (Metropolis et al., 1953; Press et al. 1992) 
is used to improve this initial TZ{m,9): 

(1) Calculate £[TZ] for this initial guess. 

(2) Make a change STZ to 71. Calculate S[TZ + 5TZ], Pi[TZ + 
5TZ] and thus £[TZ + STZ]. 

(3) If AC = C[n + STZ] - C[TZ] > accept the change STZ; 
otherwise accept it with probability exp{A£,/T), where 
the "temperature" T is defined below. 

(4) Go back to step 2. 

In this scheme each grid point {i,j) has an associated 
ATZij that governs how much TZij can change by in step 
(2). The natural initial value for all the ATZij is the RMS 
difference between Sij and the 5;^- of equation (10). In step 
(2), each change STZ is made by choosing one grid point {i,j) 
at random and adding r ATZij to TZij, where r is a random 
number uniformly distributed in (—1, 1). Whenever a change 
STZij is accepted (step 3) , ATZij is increased by a small factor 
(say 1.2); otherwise ATZij is decreased by the same factor. 
This ensures that the changes the program makes are as 
large as possible. 

Given the initial STZij, a few times rimng iterations of 
step (2) are made, and the initial temperature T is chosen 
to be the mean value of \A£\. This ensures that roughly 
half the steps will be accepted. Every Srimng iterations the 
temperature is set equal to the mean change in \A£\ of all 
the accepted steps in the previous Srimne iterations. The 
program stops with this temperature is less than ^n^ingi . 



3 KINEMATICAL EFFECTS OF KONUS 
DENSITIES 

Before investigating real galaxies, it is worthwhile pausing 
to look at the effects that konus densities can have on sim- 
ple toy galaxies. I use a Jaffe (1983) model (equation (10) 




Figure 1. The most discy (solid contours) and most boxy (dashed 
contours) density distributions p{R, z) that project to the sur- 
face brightness distribution of a flattened Jaffe model viewed at 
i = 60°. The dotted contours plot the density of the original Jaffe 
model. 

with a = 2 and P — 4) with unit scale radius, unit lumi- 
nosity and intrinsic axis ratio q — 0.6, viewed quite close to 
edge-on at i = 60°. The projected surface brightness distri- 
bution is placed on the grid (1) with n^i x ngi = 40 x 7, 
and m'l — 0.01 and m'^Q — 10. I seek the luminosity density 
p{R, z) on grid (2) with rim x ne = 60 x 10 vertices with rrii 
running from 0.01 to 100. As a test of the accuracy of the nu- 
merical projection, placing the exact p{R, z) on this grid and 
projecting yields an Sij that has an RMS difference of only 
0.05% from the results obtained by Jaffe's exact expression. 
This is much smaller than typical observational errors. 

To stop the deprojection program finding the exact so- 
lution immediately, I fix a = 1.5 in the fitting of the ini- 
tial guess (equation (10)). This initial guess has L — 0.49, 
mo = 0.26, (5 = 3.6 and an RMS error of 0.12. With 
ASij = 0.5%, Ad = 0.03, and d = -1-0.1 (discy bias) or -0.1 
(boxy bias), the program stops after about 70000 iterations 
of the Metropolis algorithm have reduced this RMS error 
to around 0.56%. To look for systematic deviations between 
Sij and Sij, I use Bender & MoUenhoff's (1987) method to 
analyze the isophotes of the solutions. This yields the semi- 
major axis and ellipticity (a, e) of the best-fit ellipse to each 
isophote, along with isophotal shape parameters 04, ae, . . ., 
obtained by expanding the angular dependence of the radial 
deviation between this best-fit ellipse and the isophote as 
a Fourier series. Both the boxy and discy solutions above 
have Ui/a ~ ±0.4%, which would be just on the threshold 
of detectability for real observations, implying that neither 
solution is acceptable. 

Letting the program continue with ASij = 0.05% in 
each case, it stops with an RMS error of 0.03% (smaller than 
the 0.05% numerical accuracy) after a further 200000 or so 
iterations. The resulting density distributions are plotted 
on Fig. 1. The projection of these densities yields isophotes 
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Figure 2. The surface brightness, ellipticity and isophotal shape 
profiles of the models in Fig. 1 viewed edge-on. I use Bender & 
MoUenhoff's (1987) definition of the isophotal shape parameters 



with \a2n/a\ < 0.03%, perfectly elliptical for all practical 
purposes. Increasing \d\ and decreasing Ad does not change 
the solutions significantly, indicating that they are the most 
extreme given the smoothness constraints. 

These densities project to the surface brightness pro- 
file of a Jaffe model for any i < 60°. So, for a fraction 
cos i — 50% of random orientations there is no sign of de- 
viations from pure ellipticity in the models' isophotes. But 
when the models are viewed edge-on, the disciness or boxi- 
ness becomes quite apparent (Fig. 2). 

The konus density obtained by subtracting one of these 
solutions from the other has almost exactly the same de- 
pendence on radius as the original Jaffe model, except at 
the very centre where a konus density can be no steeper 
than r~^ (van den Bosch 1997). Unlike the analytical konus 
densities found by Gerhard & Binney (1996), Kochanek & 
Rybicki (1996) and van den Bosch (1997), this numerical 
konus density has only approximately a cosn^ angular de- 
pendence and there are no special angles 9 along which its 
radial profile is particularly shallow. 

To investigate the kinematical effects of the photomet- 
ric degeneracy, I calculate the isotropic-rotator VPs of each 
model using the moment-based method of Magorrian & Bin- 
ney (1994). It is convenient to present the resulting VPs us- 
ing the Gauss-Hermite parameterization of van der Marel 
& Franx (1993). This gives the parameters V and a of the 
best-fitting Gaussian to each VP, along with coefficients 
hz,hi, . . . , that describe how the VP deviates from this 
Gaussian. VPs with prograde wings steeper than retrograde 
ones have h^ < 0, while VPs more triangular (flat-topped) 
than the Gaussian have /i4 > (/i4 < 0). Figure 3 shows 
the VPs of the discy and boxy models in this form, along 
with the VPs of the original Jaffe model. On the major axis 
the discy model has relatively high V and negative /13. The 



slight bias towards circular orbits that causes this also affects 
the discy model's minor-axis VPs, giving them significant 
high-velocity wings and thus raising their /14 coefficients. 
(A stronger version of this phenomenon was first noted by 
Scorza & Bender (1995) in the minor-axis VPs of the discy 
elliptical NGC 4660.) The differences between the models' 
Gauss-Hermite coefficients are approximately the same as 
the errors in high signal-to-noise measurements of real VPs. 
It is therefore plausible that one could use VPs to constrain 
the internal p{R, 2) structure of a real galaxy with a some- 
what larger disc than this simple toy model. 

4 KINEMATICAL EFFECTS OF REALISTIC 
DISCS 

To investigate the detectability of discs in real galaxies, I 
consider what the nearby edge-on E7/S0 galaxy NGC 3115 
looks like when viewed at different orientations. Scorza & 
Bender (1995, hereafter SB) have used a photometric disc- 
bulge decomposition to show that this galaxy is consistent 
with a razor-thin disc viewed at i = 84° embedded in an al- 
most perfectly elliptical bulge. The p{R, 2) profiles I obtain 
by deprojecting SB's photometry of this galaxy for i — 84° 
and i = 90° (Fig. 4) show that it can equally well be con- 
sidered as a reasonably fat disc-like component embedded 
in a slightly boxy bulge- like component. Since the galaxy is 
viewed close to edge-on, it may be possible to distinguish be- 
tween the thin- and fat-disc models using the as and higher- 
order isophote shape coefficients and a great deal of care (or, 
perhaps more sensibly, by dropping the a„ parameterization 
altogether). In what follows, let us assume that NGC 3115 
is reasonably well described by my i = 84° model, with a 
constant mass-to- light ratio and an two-integral distribution 
function. None of these assumptions are strictly true, but 
this simple model does serve to illustrate some important 
points. 

Suppose we view the galaxy close to face-on, at i = 30° . 
Fig. 5 shows that in this case the maximum value of a^/a 
is only 0.2% - there is no indication of disciness in the pho- 
tometry. The disc does, however, make itself evident in the 
model's minor-axis VPs, since the extra circular orbits seen 
from above make the VPs more centrally peaked, giving 
them large positive values of ^4 (solid curve in top panel). 
For comparison, deprojecting the i — 30° surface bright- 
ness distribution with penalty function Pshapo parameters 
d = and Ad = 4(a4/a)max ~ 0.8% yields the model with 
the most closely spheroidal isodensity contours that fits the 
data. Because of the lack of a disc in this model, it has much 
lower minor-axis h4 coefficients (dashed curve). 

The probability of observing a randomly oriented 
galaxy at inclination angle i is proportional to cosi. On the 
lower panel of Fig. 6, I plot the maximum 04/0 of the model 
of NGC 3115 as a function of cosi, showing that for about 
40% of all randomly chosen orientations (a4/a)max ^ 0.4%. 
Rix & White (1990) have found similar results using razor- 
thin exponential discs embedded in i?"^'** spheroids. As the 
disc becomes less evident in the photometry, however, its 
effect on the minor-axis VPs increases (upper panel of fig- 
ure). Indeed, there is no inclination angle for which the disc 
would not leave either a photometric or kinematic signature. 
These results are not confined to NGC 3115: I find similar 
results using SB's photometry to model the discy elliptical 
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Figure 3. The major- and minor-axis linc-of-sight velocity profiles of the models in Fig. 1 (inclination i = 60°), obtained by modelling 
each as a two-integral isotropic rotator. The panel on the left shows the Gauss-Hermite parameters of the minor-axis VPs for the 
discy (solid curves), boxy (dashed curves) and the original Jaffe (dotted curves) models. The other two panels show the Gauss-Hermite 
parameters of the major-axis VPs. 





a/arcsec 



Figure 4. (Left panel) Density distributions p{R, 2) obtained for NGC 3115 assuming i = 84° (solid contours) and i = 90° (dashed 
contours). The panel on the right shows that both models project to a reasonable fit to the observed photometry (plotted as the circles). 



galaxies NGC 3377 and NGC 4660. 



Finally, to understand how the detectability of discs 
correlates with the disc-to-bulge ratio, I embed discs with 
a range of luminosities in the toy Jaffe model of the pre- 
ceding section. Each disc has an exponential radial profile, 
with unit scale radius and an isothermal sech^(2:/2zo) ver- 
tical distribution with scale height 20 ~ 0.1. Fig. 7 shows 
that for all inclination angles the disc is clearly evident in 
either the photometry or the kinematics, provided the disc- 
to-bulge ratio is at least 3%. 



5 CONCLUSIONS 

It has been known for some time (Rix & White 1990) that it 
is possible that quite large, almost face-on discs could lurk in 
elliptical galaxies without leaving any detectable signature 
in the galaxies' photometry. More sensitive photometry is 
unlikely to improve the detectability of these discs: Gerhard 

6 Binney (1996) have shown that the deprojection of any 
axisymmetric galaxy that is not viewed exactly edge-on is 
formally degenerate to the addition of disc-like konus densi- 
ties. Using two-integral models of nearby edge-on discy ellip- 
ticals, I have shown that hidden discs do, however, leave very 
strong, easily detectable signatures in the galaxies' minor- 
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Figure 5. Evidence of the disc in NGC 3115 wiicn viewed al- 
most face-on at i = 30° . The small values of the isophotal shape 
parameters a^/a on the bottom panel show that there is no pho- 
tometric sign of the disc at this orientation. On the other hand, 
the disc leaves a strong signature in the h^ coefficients of the 
galaxy's minor-axis VPs (solid curve in top panel). For compari- 
son, the minor-axis h^ coefficients of the most closely spheroidal 
model that fits the i = 30° photometry are plotted as the dashed 
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Figure 6. Detcctability of the disc in NGC 3115 as a function 
of inclination angle i. The maximum values of ai/a for each in- 
clination are plotted on the bottom panel, showing how quickly 
the disc becomes photometrically invisible with decreasing incli- 
nation. The top panel shows the effect that the disc has on the 
galaxy's maximum minor-axis ^4 coefficient (solid curve) for each 
inclination. For comparison, the maximum minor-axis /i4 coeffi- 
cients of the most closely spheroidal two-integral model that fits 
the observed photometry are plotted as the dashed curve. 

axis VPs, provided the disc-to-bulge ratio is greater than 
about 3%. 

These results are based on two-integral models, but they 
have important implications for more realistic three-integral 



Figure 7. Detectability of the disc in the toy galaxy model de- 
scribed in the text for disc-to-bulgc ratios of 0.1 (solid curves), 
0.03 (dashed curves) and 0.01 (dotted curves). As in Fig. 6, the 
bottom panel shows the maximum 04/0 coefficient as a function 
of inclination, while the top shows the maximum minor-axis hi 
coefficient. 

models. The most sophisticated current modelling machin- 
ery (e.g., van der Marel et al. 1998; Gebhardt et al. 1998) 
uses an assumed (usually spheroidal) p{R, z) and corre- 
sponding self-consistent potential to constrain a galaxy's 
orbit distribution from its VPs. Given, say, an isotropic 
galaxy with a weak, invisible face-on disc, it is possible that 
such machinery could fit the positive hi coefficients of the 
galaxy's minor-axis VPs using a spheroidal p{R, z) with a 
radially biased orbit distribution. (In principle, a particu- 
larly cold disc might leave an unambiguous signature on the 
VPs, but it is unlikely that this signature would survive 
the instrumental broadening that real measurements suffer 
from.) Thus it is important that some attention be focused 
on modelling galaxies that are known to be close to edge-on 
(i.e., are discy), so that the uncertainty in the photometric 
deprojection is minimized. 

Ultimately, however, one is interested in all the other 
galaxies whose inclination angles, density distributions and 
orbit distributions are unknown. Until more general mod- 
elling machinery that "knows" about the orbital structure of 
edge-on galaxies becomes available, it might be worthwhile 
relaxing the self-consistency requirement in current models 
and have them fit to VPs and photometry, rather than fitting 
to an assumed, unobservable p{R, z). It would be possible to 
iterate, using the potentials of the resulting density distribu- 
tions to generate proper self-consistent models. This could 
be an informative first step in trying to constrain both the 
internal shapes and orbital structures of galaxies. 

Unfortunately, it is not possible to test these ideas yet 
because there are currently very few power-law galaxies with 
high signal-to-noise measurements of minor-axis VPs. Let us 
hope that this is only a temporary problem. 
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sight, I choose n — 4na. 

Since p{m, 6) is obtained by linear interpolation in 
TZ{m,6) = log p{a, 6) this sum can be rewritten in the form 



, Teller A., Teller 

Vctterling W.T., 
Cambridge Univ. 



127, 
Dor- 



(15) 



} Ek exp (^AkTZak,bk + BkTZak + l,b^ 
fc=0 

+ CfcTlafc.bfc + l + -Dfc7^afc + l,6fc+l) , 

where the (a;',2y')-dependent coefficients Ak-.-E^, ak and 
bk do not depend on 7?. and so only have to be evaluated 
once for each (a;^^-, j/^^). Expressions for these constants are 
quite easy to work out, but tedious to write out. 

If, like Romanowsky & Kochanek (1997), I had cho- 
sen to interpolate linearly in p rather than log p, then equa- 
tion (15) would be replaced by an even simpler expression. 
However, this interpolation-in-p scheme is unable to repre- 
sent steeply cusped density distributions accurately: using 
it, the RMS error in the numerical projection of the toy 
model used in Section 3 is 1.6%, rather than the 0.05% ob- 
tained with interpolation in log p. On doubling Um to in- 
crease the resolution of the radial rrii grid, the error goes 
down only slightly to 1.3%. (On the other hand, the error 
in the logp scheme goes down to 0.015%.) Thus, the extra 
expense that interpolating in logp brings to the numerical 
projection equation (15) is easily justifiable. 
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APPENDIX A: THE PROJECTION ALONG 
THE LINE OF SIGHT 

The projected surface brightness distribution I{x' , y') of the 
density distribution p{R, z) is 

/oo 
p(R,z)Az' (11) 

-oc 

where {R, z) is related to {x\y',z') through the rotation 

[ y J = [ cosi sinj ] I v' ] , (12) 

\z / \ sin i — cos i / \z' / 

and E? — x^ + y'^. With the substitution 

/2 , /2\l/2 



z = (x + y ] smhu, 



(13) 



the integral in equation (11) can be approximated using the 
trapezium rule as 



I{x', y') ~ {x''^ + y'^f'^ du ^ p{Rk, Zk) cosh 



Mfc, 



(14) 



with U{_) and u„ chosen such that the corresponding z' lie 
just outside the edges of the model density grid, and the 
intermediate Uk spaced linearly between them with du = 
(u„ — uo)/n. The points Rk and Zk are the intrinsic co- 
ordinates of the point {x' , y' ,z'k) in the fc*'' term in the sum. 
To ensure fair sampling of the grid points along the line of 



